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METHOD AND APPARATUS OF RELATIVE DATAPATH CELL PLACEMENT 
WITH STRUCTURE BONDING AND EXPONENTIAL FIELD SAMPLING 



The invention relates to computer-aided design (CAD) of integrated circuits, and 
particularly to placement of cells in datapath design. 

Description of Related Art 

T^wo general categories of placement methods are used, the first one random 
10 placementVid the second one is structural placement. For random placement, cells are 
connected byUist of nets. The objective of placement is to minimize a cost function that 
are related to total sum of net length with fixed constraints of placement area, routability, 
and timing. sWeral conventional placement algorithms relating to this type of 
placement method hVe been disclosed. 
15 In bonding cell instances, a conventional solution places wires at locations 

between cell instances which produce shortest net lengths. Such approach produces an 
undesirable random or irregular structure arisen from the placement of uncorrelated cell 
instances. Another conventional solution places cell instances in a rigid structure by 
matching pin locations between cell instances. This technique is prohibitively restrictive 
20 since adjacent cell instances must be aligned with matching pin locations and hence 
introduces many dead placement spaces. 

In a density map partitioning, conventional solution adopts a uniform grid 
partitioning. This approach requires substantial computation time in calculating the 
density forces of a given point on a partition map. Because each square or rectangle in 
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the grid is uniform in size, the computation time is greater if the grid size of a cell 
instance is smaller which produces a high number of cells. 

Accordingly, it is desirable to place cell instances in an optimized structure 
bonding. It is further desirable to have a density map partition which is non-uniform in 
computing a force update vector in reducing the extensive computational complexity 
time. 





SUMMARY OF THE INVENTION 

The invention discloses a relative placement of datapath elements by keeping the 
relativitiyty of ce\l instances in a datapath group without rigid alignments between cells. 
10 In one embodiment to encourage placement of a desirable structure, pseudo cells, pseudo 
pins, and pseudo nets are selected to be placed at certain locations with respect to real cell 
instances. The end result produces a cluster of real cell instances that form a desirable 
structure while minirmzing the length of nets. The invention therefore achieves 
objectives of minimizinghhe sum of net lengths while encouraging cells instances to form 
15 in a desirable structure, sush as a column structure, a row structure, or a square structure. 

In a fuither aspect of the invention, a non-uniform partitioning of a density map 
for calculating amorce update vector is disclosed. The partitioning is taking over a region 
A to compute Rieirann sum approximations of a function F over the region A. A force 
update vector is calculated for a given cell instance within the region A where 
20 neighboring cell instance&have an exponentially larger grid size as cell instances extend 
further away from the givenxell instance. 
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Advantageously, the structure bonding reduces the dimension of a chip by 
minimizing dead placement spaces in a datapath structure for efficient utilization. 
Moreover, the structure bonding minimizes data skew of datapath. The structure bonding 
of cell instances with pseudo elements also encourages low timing delays. Additionally, 
the exponential portioning advantageously reduces significantly on the computational 
time in the determination of a force update vector on a given cell instances arisen from 
attractive and repelling forces. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a flow diagram illustrating the general process in the placement of cell 
instances in accordance with the present invention. 

FIGS. 2A-B are pictorial diagrams illustrating relative cell placement of a column 
structure and a row structure in accordance with the present invention. 

FIG. 3 is a pictorial diagram illustrating relative datapath placement in accordance 
with the present invention. 

FIG. 4A is a pictorial diagram illustrating cell instances constructed with a pseudo 
net and pseudo pins in accordance with the present invention; FIG. 4B is a graphical 
diagram illustrating cell instances with relative positioning in accordance with the present 
invention. 

FIG. 5 is a pictorial diagram of one embodiment in field driven placement with 
positive and negative charge particles in accordance with the present invention. 

FIG. 6 is a flow diagram illustrating the process in a structure bonding in 
accordance with the present invention. 

-4- 

Structure and Sampling ARC-P002 



# 



*! 3 



■Q 10 
IB 



□ 15 



FIG. 7 is a flow diagram illustrating the process of a Riemann sum approximation 
of the force update vector in accordance with the present invention. 

FIG. 8 is a pictorial diagram illustrating one embodiment of an exponential field 
sampling partition in accordance with the present invention, 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENT 

FIG. 1 is a flow diagram of illustrating process steps by a placement engine 10 for 
the placement of cell instances in a structure bonding. Placement engine 10 processes 
placement steps so that each step generates a new placement based on an existing 
placement. Each placement step modifies an existing placement by combining attraction 
actions between connected cell instances with repulsive actions between cell instances. 
Placement engine 10 receives 11 an input netlist and identifies 12 cell structures from the 
netlist. Cell structures are extracted 13 from the netlist by placement engine 10 to 
perform 14 bus analysis. Placement engine 10 then extracts 15 general regularity 
structures and creates 16 bonding between the extracted cells. Subsequently, placement 
engine 10 performs 17 of field-driven placement of cell instances on a density partition 
map. 

FIG. 2A is a pictorial diagram illustrating one fundamental relative structure in 
dapapath in a column structure 20. Column structure 20 embodies a vertical sequence of 
cell instances being fixed while horizontal distance between two cell instances is not 
fixed. Preferably in a column structure 20, cell instances are within the vicinity with one 
another. For example, in an 8-bit datapath, the binary "0" bit cell is placed at the bottom. 
The binary "1" cell is then placed above the binary "0" bit cell. This vertical sequence 
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continues with the binary "2" bit cell being placed on top of the binary "1" cell until all 

eight cells are stacked in a chronological sequence from the lowest binary value to the 

highest binary value. Although the 8-bit vertical sequence in the datapath is fixed, cells 

that are horizontally adjacent to the 8-bit datapath have the flexibility to be placed in the 

vicinity, without being strictly aligned, to the eight vertical stacked cells. FIG. 2B is a 

pictorial diagram illustrating another fundamental relative structure in dapapath in a row 

structure 25. Row structure 25 embodies horizontal sequence of cell instances being 

fixed while vertical distance between two cell instances is not fixed. A further 

fundamental relative datapath structure is an array structure where both the horizontal 

sequence and vertical sequence are fixed. One skilled in the art should recognize that 

other type of structures or modifications can be practiced without departing from the 

scope of the invention. 

Some terminologies used in describing the invention are defined as follows. A 

cell instance is an object which inhibits, typically a rectangular subset of a placement 

area. A net is an object which connects some set of cell instances. A pin is an object 

which represents the relative location in a cell instance's region where a net connection 

occurs. A placement is a mapping from a set of cell instances to a set of positions in a 

placement area. These positions can represent, for example, the centers or lower left 

corners of regions that cell instances inhibit. A placement area represents some subset of 

a two-dimensional real space, which is typically a bounded rectangle or rectilinear region. 

A total wire length in a placement is the value of some function that approximates the 

cost associated with the physical wire that connects to physical objects in which cell 

instances represent. The write length function is chosen to minimize the sum of half- 
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perimeters of net bounding boxes. For each net, a minimal rectangle is found with 
horizontal and vertical edges that contain all pins associated that net, and add the width 
and height of that rectangle to a total wire length function. The utilization ^ of a 
placement is the sum of areas of cell instances, divided by an area of the placement area. 
The utilization is typically between zero and one, and is independent of the positions of 
cell instances. A value of zero means that there are no cell instances, and that the 
placement area is empty. A value of one means that there is exactly enough room for all 
cell instances, and there would be no unused space if a cell instance can be placed without 
overlap. 

FIG. 3 is a pictorial diagram illustrating a structure bonding 30 of cell instances. 
Structure placement is to place a particular group of cell instances at certain locations that 
are relative to one another. Cell instances 31a-g on the left side are concatenated to cell 
instances 35a-g on the right side in a rectangular structure, even though cell instances 
31a-g are not perfectly aligned with cell instances 35a-g. Placement engine 10 places 
cells instances 31a-g and 35a-g with the aid of pseudo nets thereby producing structure 
bonding 30. Structure bonding 30 results in a zig-zag structure, the usage of pseudo nets 
encourage cell instances to form a structure without introducing extra dead placement 
spaces. 

FIG. 4A is a pictorial diagram illustrating cell instances constructed with pseudo 
net and pseudo pins. The shortest distance between a cell instance 40 and a cell instance 
41 is a straight line from the top of cell instance 41 to the bottom of cell instance 40. 
Therefore, a pseudo pin 42 is placed in bottom of cell instance 40 for connecting with a 
pseudo pin 43 at the top in cell instance 41. A pseudo net 44 connects between pseudo 
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pin 42 and pseudo pin 43 together while minimizing the length of wire interconnect. 
Although pseudo pins 42, 43 and pseudo 44 do not actually exist, pseudo net 44 and 
pseudo pins 42, and 43 are treated as real nets and pins for encouraging the placement of 
cell instances in a structure. 

A pseudo cell instance c f 48 is added with a width and height equal to zero, 

without interfering with a density map. The pseudo cell instance c f 48 is assigned to a 

position on the density map to encourage a desirable structure. Between pseudo cell 

instance, c f and real cell instances c„ c 2 , ... c^, a set of d\ 9 dk relative positions is 

created and satisfies: xj - x/ + dj for all j, where the term x/ denotes an optional 

pseudo instance and the term dj denotes a relative position. 

An example is shown in FIG. 4B. A set of thee cell instances c t 45, c 2 46, and c 3 
47, each represents a square of five units in length and width. The three cell instances c, 
45, c 2 46, and c 3 47, are aligned in a cascaded row. If the lower left corner of c x 45 is at a 
location jti, then cell instance c 2 46 has a lower left corner at x 2 = x\ +(5,0) , and cell 
instance c 3 47 has a lower left corner x 2 =x\ +(10,0). In this example, pseudo cell 
instance c f 48 is created for sharing the location x\ with the lower left corner of c,. 
Induce relative position are computes as d\ = (0 ,0) between c f 48 and c, 45, d 2 = (0 ,5) 

between c f 48 and c 2 46, and di = (0 , 10) between c f 48 and c 3 47. 

Upon determining the relative positions between the new pseudo cell instance and 

real cell instances, new pseudo nets are created to bond the new pseudo cell instance and 

real cell instances. Although several techniques are used to construct pseudo nets, the 

-8- 

Structure and Sampling ARC-P002 



preferred technique is to select a connection which minimizes wire lengths. For example, 
a simple net is used to construct cell instances in a structure. If, however, large nets are 
given lower weights than small nets, then the bonding effect gained by the net may be 
diminished. In such scenario, two-pin nets may be created among cell instances. The 
pins of nets are placed at locations that encourages a desire relative placement of the 
structure. 

In assigning relative locations of pins to their corresponding cell instances, pins 
associated with the pseudo or fake cell instance c f 48 are placed at a relative location of 

(0, 0) from the position of c f 48. Subsequently, a pin that is associated with a real 
instance Cj may be assigned to a relative position of- d, from c ; . If a relative location of 
(0, 0) is not assigned to the pseudo cell instance's pin, then the difference of these relative 
locations is d p i.e. p f - p 7 = d y where p f represents the location ofc / 's pin relative to 

x p while pj represents the location of c } 's pin relative to x y . 

Continuing with the example above, three fake nets n l5 n 2 , and n 3 which connect 
c f to , respectively, c ; 45 c 2 46 and c 3 47. If a pin is placed for each net at a location of 

(0, 0) relative to the pseudo cell instance, a pine for n, is place at -d x = (0, 0) relative to 
c 1? a pine for n 2 is place at - d 2 = (0, -5) relative to c 2 , and a pine for n 3 is place at -d 3 = 
(0,-10) relative to c 3 . 

Preferably, placement engine 10 attempts to make pseudo nets as short as possible 
between cell instances. The shortest possible net is one with a length of zero, in which all 

of the net's pins are at the same position z . For example, if Cj is a real member of such 
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a net, and if pins are placed at p, relative to c y and p f relative to then a desirable 

result is achieved x y = Xy+ p r - p, = Xy+ d y , where z= x y + p y = x 7 + p r . The 
optimal structure is to have all pseudo nets have zero length so that cell instances in a 
structure have exactly the desirable relative positions. 

Optionally, rather than constructing pseudo cell instances, pin offset vectors can 
be assigned as if a pseudo instance existed. Each cell instance in the structure is 
connected to surrounding cell instances by a path of pseudo nets. 

Members of a structure may be scattered if real nets connecting to cell instances 
overpowers pseudo nets. In such scenario, placement engine 10 increases the weight of 
pseudo nets to reduce the effect of undesirable structure. In one embodiment, the weights 
of pseudo net are set equal to the weights of real nets to maintain structure cohesion. 

FIG. 5 is a flow diagram illustrating the process of structure bonding 50 in 
accordance with the present invention. A pseudo instance is created 51 and selects 52 at a 
relative location to encourage structural placement. Placement engine 10 finds 53 offsets 
between pseudo instance and real instances, represented by the symbol d y 
Subsequently, placement engine 10 selects offset between pseudo instance and the 
associated pseudo pin, p f . After pin offsets are selected 55 for real instances and their 
associated pseudo pins, p j9 placement engine 10 creates a pseudo net connecting pseudo 
pins. 

FIG. 6 is a pictorial diagram of one embodiment in a field driven placement 60 

with positive and negative charge particles. Particles 61 and 62 are tied together with a 

spring 63. Similarly, particles 64 and 65 are connected with a spring 66. Also, a spring 
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69 ties particles 67 and 68. Placement engine 10 balances different objectives through a 
sequence of placement steps. Each placement step modifies the placement by combining 
attractive actions between connected cells with repulsive actions between cell instances. 
Attractive actions serve to optimize objective functions while repulsive actions serve to 
reduce overlap among cell instances. The name "field driven placement" encompasses 
positively charged particles placed in a negatively charged field. The particles are tie 
together with springs. The objective is to minimize the following function: 

min YsW 2 + -jE 2 dr Eq. (1) 

The term r 2 in Equation 1 represents the energy of a spring, which can also be 
written as (x { - x } f + (y { - y^) 2 . The additional term of ~ ^E 2 dr represents the energy of 

the field. To minimize Equation 1, a vector differentiation is taking on the expression on 
all movable objects location. The derived Equation 1 becomes: Ax + 6+ F x = 0, and 

Ay + 6+ = 0, where (Ax + b, Ay + b) represents the spring force and (F x , F y ) 

represents the electrical force. 

The derived Equation 1 can be interpreted in which the sum of an overall force is 
equivalent to zero at equilibrium stage. The electrical force is then expressed as follows: 

Wx>y>^(x>y)+K(x>y)™& fyx>y>^{x>y)+K( x >y)- The term f»(x, y) 

represents the force due to attractive force of the positive plane, while F xr represents the 
repelling force from other cell instances, and thus 
Wx,y) = ^(x 9 y) 9 F^(x 9 y) 9 y r (x 9 y) = {F xr (x>y\F yr (x,y)} . Cell instances move in a 
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plane. According to Gauss's Law: e dEdS -q^> E = 



. The attraction force F, 



a 



2/rer 




the left and right boundaries in a placeable area, and b and t represents the top and bottom 
boundaries in the placeable area. Similarly, the repulsive force acting upon a cell is 
5 computed as the summation of the repulsive force from other cells: F r (x,y)& 



An force update vector is an accumulation of repulsive actions upon a given cell 
instance by other cell instances. These actions are determined such that cell instances are 
pushing away from high density regions while pulling towards low density regions. The 



U 10 constant value of k affects the strength of the force update vector. Low values of k 

iii 

W reduce the tendency of placement engine 10 to favor overlap reduction over the 

a optimization of objective functions, with a tradeoff against runtime in computing a 

greater number of force update vectors. High values of k presents a converse issue, in 

which a poor solution is achieved quickly. 
15 Another way to characterize the placement step and the computation of an force 



update vector is described as follows. An initial placement step is to update a force 
vector / c for each cell instance c, which is used as an offset vector when finding that 
instance's location. The force vector / c , which is initially set to zero, is additively 
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updated during each placement step by repulsive actions from other cell instances, as well 
as attractive actions from the placement area. 

The subsequent placement step comprises an iterative process of a small number 
of intermediate updates of placement. Each intermediate updates include two phases, one 
5 which determines virtual positions of the nets and one which determines the positions of 
the cell instances. The first phase of the intermediate update finds a virtual positions x n 
for each net n which connects a set of cell instances {c l5 . . Cj). The position is obtained 
by computing a weighted average of instances' positions: 

10 where {w„ Wj) are non-negative cell instance weights. A constant weight of 1 is 
typically used. 

During the second phase, an intermediate position x c is found for each cell 
instance c. This position is a shifted weighted average of the virtual positions of the nets 
{n 1? . . n k } which are connected to c. The weighted average is shifted by the force vector 
15 which is computed at the outset of placement steps: 

- l7 

x o + A 

2 k 

where {w l5 w k } are nonnegative net weights. Typically, the weight w^s some (— ) 

P 

where p is the pin count of the net n^ (i.e. the number of cell instances which are 
connected to nj, and k is some non-negative integer. By starting with k = 0 and 
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gradually increasing to k = 2, the total overall wire length is minimized when making 
transition from a randomly chosen placement to another placement. 

Optionally, each net is associated with a force vector which is used as an offset for 
that net's virtual position, and cell instance positions receive no such offset. This 
5 method is referred to as dual placement. 

At each placement step, the force vector associated with each cell is updated. The 
update involves a multiplication by a force trimming scalar first, followed by an addition 
with a force update vector. The force trimming scalar is typically some real number 
□ between zero and one and determines the balance between competing goals of 

! P 10 discouraging overlap and optimizing objective functions. Values that are close to 0 
encourage shortening nets, while values that are close to 1 encourage cell instances to 
;S spread out. The force update vector represents an accumulation of repulsive actions upon 

O a given cell instance by other cell instances. Repulsive actions are determined so that cell 

|1J instances are pushed away from high density regions and pulled towards low density 

G 15 regions. The force update vector is represented as follows: 



where k represents some constant, and D represents a density function. The 
integral is taking over an entire placement area, or all of two-dimensional space as long as 
the density function is uniformly zero outside the placement area. The force update 
20 vector e c associated with a cell instance c is then given by e(x c ), where x c is the 
position of c. The constant value k affects the strength of the force update vector. A low 
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value of k reduces the tendency of placement engine 10 to favor overlap reduction over 
optimizing objective functions at the expense of runtime. A greater number of force 
update vectors is computed. A high value of k have the opposite concerns, which 

3 

achieves poor solutions quickly. In one embodiment, a value of is used where (i 

1000-// 

5 represents the utilization of the placement. 

A density function D maps a given position x ? in the placement area to a value 
which is given by the number of cell instances which cover the position x\ less the 
utilization of the placement: D (x ) = ^Z c ( x ) where ^represents the characteristic 

c 

function associated with the region that instance c inhabits, i.e. Z c ( x ) equals one if the 

10 position x ' is covered by c, and otherwise equals to zero. The subtraction of the 
utilization from the sum of the characteristic functions serves to provide an attractive 
action upon cell instances towards regions with low density. 

The integral in the force update vector is computed through a Riemann sum 
approximation. FIG. 7 is a flow diagram of this process 70. The placement engine 10 

15 finds 71 the positions of the sample regions and selects 72 a region. Within that region, 
the placement engine 10 selects 73 a given point to serve as a representative for 
computing the part of the force update vector which represents attractive and repulsive 
actions in that region. The placement engine 10 adds a scaled effect of that point to the 
overall effect. If there are more regions, the process 70 returns to the operation 72. 

20 Otherwise, the placement engine 10 completes the process 70 in a step 76. 
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To approximate the integral of a function F over some region A, the region A may 



be partitioned into a collection of pairwise disjoint subregions A l9 A k , and then 



representative points & x 



, . . ., a k may be chosen so that aj € A } for all j = 1, . . ., k. If |4/|| 



represents the area of region Aj , then the corresponding Riemann sum approximation 



5 



{/is: t^jAjl 



Riemann sums are typically computed by dividing a region A into a rectangular 



grid, where subregions occupy an equal area. Representative points can be chosen to be 



comers, centers, random locations, or places where the function to be integrated is 



— 



minimized or maximized. To make the computation of the force update vector more 



W 10 efficient, the placement area is divided into such a uniform grid, and store a density map 
^ y that indicates, for each uniform region in the grid, the cell density within that region. 

This partition is called the density map partition. 

11! 

Ifi There are some drawbacks to using a uniform grid of regions in the computation 

O 

O of the force update vector. Selecting a small grid size promotes accuracy in Riemann sum 

15 approximation, but with a severe cost in computation time. Conversely, selecting a large 
grid size reduces the computational time of force update vectors, but allows for a large 
amount of irregularity in the cell density within the regions. For example, if 10,000 cell 
instances are to be placed in an area that is partitioned into two subregions, with 5,000 
cell instances for each subregion, then the placement may be categorized as uniform if 
20 each of the regions has its 5,000 cell instances occupying the representative position 
chosen for that region. 
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To avoid the difficulties arising from the use of a uniform partition in the 
Riemann sum approximation, a non-uniform partition is chosen whose members are 
made up of groups of smaller density map grid elements. To this end, a convenient 
labeling system is established for the density map partition elements base on rows and 
columns. Given a model rectangle E, the two-dimensional plane is tiled with a grid of 
copies {Ey I yeZ} of E, with the following order: i 0 < i x if and only if E kJ is left of 

E hJ for all j; and j 0 < j, if and only if E.. if below of E iJx for all i. Without loss of 

generality, the E 0 0 region may be placed at any location on the placement area, such as at 

the lower left corner. 

The density map partition induces a map 8: Z 2 -> 91 such that S(i, j) is equal to 
the cell instance density in E itj . For example, if E 4 7 contains 3 area units of cell instance 

c 1? and 2.4 area units of cell instance c 2 , and no other cell instances intersect E 4 7 , then 
5(4, 7) = 3 + 2.4 = 5.4. This map, less the utilization (i, is used as the density function D 
in the Riemann sum calculation of the force update vector. 

In describing a non-uniform partition used in Riemann sum approximation, two 
families of intermediate partitions are described to induce a final partition by way of 
intersections of intermediate partition members. The intermediate partitions divide a 
plane into either horizontal or vertical stripes, with a central thin strip that is surrounded 
by increasingly thick stripes. 

Given a density map partition {E & \ ijeZ}, a base index i 0 e Z , and m € Z , 

V io (m) can be defined as, the m* vertical stripe about i 0 , as follows: 
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V i AO)=u{E io JjeZ} 

V. (m)= yj{E. t J e Z, and 2 m_1 <k < 2 m }for m>0 
V, (-m) = u{£. ,.jeZ, and 2 m_1 <k < 2 m }for m>0 

Similarly given a base index j 0 eZ, and n gZ, H h (n) can be defined, the n* 

5 horizontal stripe about j 0 , as follows: 
H jo (0)«u{£,J/eZ> 

H h («)= u {E ijo+k \j e Z, and 2"" 1 <k < 2" } for n>0 
H jo (-«) = u {E i jg _ k | i € Z, a«J 2"' 1 <£ < 2" }for n > 0 

To verify the families of vertical stripes and horizontal strips: let i 0 e Z. The 



~ : 10 collection { V { (m)\ m € Z) is a partition of a plane. For m > 0, the width of V i (m + 1) is 

□ 

jf; twice that of V k (m) , and 2 m times that of (0) . For m < 0, the width of V k (m - 1) is 

□ twice that of V k (m) , and 2" m times that of V io (0) . 

In the Riemann sum calculation of a force update vector, the exponential partition 
is defined as follows. Given a base index 0 0 ,y 0 )eZ 2 , and integers m and n, 
15 A iojo (m,n) can be defined, the (m, n) ^ exponential partition member, as 
V io (m)nHj Q (n). If there is no uncertainty about the values of i 0 and j 0 , then the 

abbreviation of A (m, n) is used. 

The exponential partition members are so named because their areas increase 

exponentially as the indices m and n move away from zero. Let (i 09 j 0 ) e Z 2 . Then, 
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A io jo (0, 0) = E iQ fjo . The collection {A (rn, n) \ 



m,neZ} is a partition of the plane. For m 



h 10 



> 0 and n > 0, the area of A(m + 1, n + 1) is 2 m+n times that of A(0, 0). FIG. 8 is a pictorial 
diagram of one embodiment in implementing an exponential field sampling partition. 

The force update vector of a given cell instance c is then computed as described 
below. Let E iQ9j represent the density map partition that contains the center x of c. 

Representative locations are selected { x m J m 5 n € Z } so that x m n e A i0 , j0 (m, n) for all 
m and n in Z. For each m and n in Z, the integers i m and j n are found such that x m n 
€ E. . . A term D(m, n) is defined to be zero for m and n in Z if E- . does not 



intersect the placement area; otherwise, D(m, n) is define as ^(i m? jj - //, where 
^represents the density map and //represents the utilization of the placement. The 
resulting Riemann sum approximation is as follows: 



Preferably, representative locations are selected randomly to avoid undesirable clustering 
of the cell instances. 

FIG. 9 is a pictorial diagram of an example of one embodiment in implementing 
an exponential field sampling partition. A given cell instance, denoted by the symbol 
is selected, along with its region in the density map which will serve as a central region. 
Each of the eight regions immediately surrounding the point # is a lxl region unit. 
Moving outward, the size of a region increases exponentially, from 2 units, to 4 units, and 
more. The farther away a square or rectangle is from the point *, the larger the sampling 
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area. Because a single sampling point is chosen for each region, the larger (and thus 
farther) regions take less processing time per unit of area than the smaller (and closer) 
regions do. The choice of a single sampling point for a larger region, for example a 4x4 
region, also reduces the accuracy of the estimation of that region's effect in the force 
5 update vector, but this lack of accuracy is offset by the fact that such regions are far from 
& and, because of the inverse-distance-squared term in the force update vector, have a 
smaller effect per unit of area on * than a closer region, such as a 2x1 region. 

The above embodiments are only illustrative of the principles of this invention 

□ and are not intended to limit the invention to the particular embodiments described. For 
=P 10 example, variations of pseudo elements may be used to encourage cell instances to be 
I*"; placed in a structure bonding. Furthermore, the exponential partitioning of a density map 
[% may have various assigned ratios that are not necessarily powers-of-two. Accordingly, 
p various modifications, adaptations, and combinations of various features of the described 
j tj embodiments can be practiced without departing from the scope of the invention as set 

□ 1 5 forth in the appended claims. 
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